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ABSTRACT 

Context. A study of the structural and scaling properties of the temperature distribution of the hot. X-ray emitting intra-cluster medium of 
galaxy clusters, and its dependence on dynamical state, can give insights into the physical processes governing the formation and evolution of 
structure. 

' Aims. Accurate temperature measurements are a pre-requisite for a precise knowledge of the thermodynamic properties of the intra-cluster 
,: medium. 

^\ Methods. We analyse the X-ray temperature profiles from XMM-Newton observations of 15 nearby {z < 0.2) clusters, drawn from a statistically 
representative sample. The clusters cover a temperature range from 2.5 keV to 8.5 keV, and present a variety of X-ray morphologies. We 
• derive accurate projected temperature profiles to ~ 0.5i?200. and compare structural properties (outer slope, presence of cooling core) with a 
' quantitative measure of the X-ray morphology as expressed by power ratios. We also compare the results to recent cosmological numerical 
, simulations. 

rS Results. Once the temperature profiles are scaled by an average cluster temperature (excluding the central region) and the estimated virial 
radius, the profiles generally decline in the region 0. 1 R20Q < R < 0.5 Rioo- The central regions show the largest scatter, attributable mostly to the 
presence of cool core clusters. There is good agreement with numerical simulations outside the core regions. We find no obvious correlations 
between power ratio and outer profile slope. There may however be a weak trend with the existence of a cool core, in the sense that clusters 
^ ' with a central temperature decrement appear to be slightly more regular. 

• • I Conclusions. The present results lend further evidence to indicate that clusters are a regular population, at least outside the core region. 
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1 . Introduction profiles are typically determined with considerably less spatial 

resolution than density profiles. 

The temperature and density are the key measurable charac- ^j^^ measurement of radial temperature profiles is further 
teristics of the hot. X-ray emitting intracluster medium (ICM). complicated by the density squared (n^) dependence of the X- 
The determination of important derived properties such as en- emission. The steep drop of the X-ray surface brightness 
tropy, pressure, and, under the assumption of hydrostatic equi- ^^^j^ ^-^^^^^^ ^^^^ ^j^^ ^^^^^^^ combined with the background 
librium, the total mass, is dependent on accurate estimation of ^^^^ ^^^^-^^ ^^j^^ instrumental sources, makes accurate 
these quantities. Because of limited photon statistics' it is usual ^leasurement of the temperature distribution at large distances 
to measure the density and temperature in terms of radial pro- ^^^^ ^^^^^ ^ technically challenging task, 
files. However, while the density of the ICM is relatively easy 

c ^. f u • u» ci f ■ 1 The earliest temperature profiles were mea- 

to measure from the surface brightness profile of a given cms- ^ ^ 

^ » J » • ^- ■ fc ■ . u » sured with Einstein, EXOSAT, Spacelab-2 and 

ter, the temperature deternunation requires sufncient photon ' 
. . u ij J . T-u Tr^\/i t » GINGA only for the nearest, brightest clusters (e.g., 

statistics to build, and fit, a spectrum. Thus ICM temperature ^ , . t n ^ • , , • n ^ . 

Fabricant, Lecar & Gorenstein 1980; Fabricant & Gorenstein 

1983; Hughes, Gorenstein & Fabricant 1988; Eyles et al. 
5'endo^n«/regMeif^to:G. W.Pratt, e-mail: gwp@mpe.mpg.de 1991; Koyama, Takano & Tawara 1991). The low, stable 

' Also the need for an azimuthally symmetric approximation for background of ROSAT made possible spatially resolved 
purposes of deprojection. spectroscopy of poor clusters (e.g. David, Jones & Forman, 



2 



G.W. Pratt et al.: X-ray cluster temperature profiles 



1995); however, limited spectral resolution and band- 
width made such measurements difficult for hotter clusters 
(e.g. Henry, Briel (& Nulsen, 1993; Briel & Henry, 1994; 
Henry & Briel, 1995). ASCA and BeppoSAX had sufficient 
high-energy sensitivity to accurately measure the temperatures 
of hot clusters. However both of these satellites suffered from 
significant PSF blurring, which, in the case of ASCA, was 
exacerbated by a significant energy dependence. As a result, 
at the end of the ASCA/BeppoSAX era, the exact shape of 
cluster temperature profiles was still under vigorous debate 
(Markevitch et al., 1998; Irwin, Bregman & Evrard, 1999; 
White, 2000; Irwin & Bregman, 2000; Finoguenov et al., 
2001; De Grandi & Molendi, 2002). 

Chandra and XMM-Newton do not suffer from major PSF 
problems. The on-axis Chandra PSF is negligible, while the 
XMM-Newton PSF becomes an issue only for clusters with 
very centrally peaked core emission; in addition, neither is 
energy-dependent. Recent observations of moderately large 
samples consisting primarily of nearby cooling core clus- 
ters wiffi XMM-Newton (Piffaretti et al., 2005) and Chandra 
(Vikhlinin et al., 2005) have largely validated the original 
ASCA results of Markevitch et al., which suggested that tem- 
perature profiles declined from the centre to the outer regions. 
However, other Chandra and XMM-Newton observations have 
found flatter profiles (Allen et al., 2001; Kaastra et al., 2004; 
Arnaud et al., 2005). As of the time of writing, no systematic 
attempt has been made, with either XMM-Newton or Chandra, 
to look at the temperature profiles of a representative sample 
of nearby clusters^. Although other projects on representative 
samples are in progress (e.g., Reiprich et al., 2006), they are 
not expected to be able to map the temperature distribution out 
to large radius. 

In this paper we deal with observations of 15 clus- 
ters from a statistically representative sample observed with 
XMM-Newton . We describe in detail the data reduction and 
background subtraction, and compare our results with previous 
work and with those from cosmological hydrodynamical sim- 
ulations. We also make a preliminary investigation of corre- 
lations with quantitative morphological measures. We present 
only projected temperature profiles in this paper - such pro- 
files are direct observables and do not depend on complicated 
PSF and deprojection algorithms. We will deal with correc- 
tion of the profiles in forthcoming papers which make use of 
observations of the full sample. All results are given assum- 
ing a ACDM cosmology with Q,„ = 0.3 and Qa - 0.7 and 
Hq - 70 km s"' Mpc"'. Unless otherwise stated, errors are 
given at the 68 per cent confidence level. 

2. The sample 

The XMM-Newton Legacy Project for the study of cluster 
structure was initiated to study the structural and scaling prop- 
erties of a large, representative sample of clusters. Since full 
details will appear in a forthcoming paper, we present here only 
a short summary of the sample selection. 

^ Some work has been done on medium-distant clusters, see 
(Zhang et al., 2004; Kotov & Vikhlinin, 2006) 
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Fig. 1. X-ray luminosity - redshift ( Lx - z) distribution of the 
REFLEX cluster sample in the redshift range < z < 0.25. The red- 
shift/luminosity selection for the Legacy Project sample is indicated 
by the boxes. Filled circles indicate clusters from the Legacy Project 
which are discussed in this paper; dotted circles indicate clusters for 
which re-observations are necessary. The solid line is the REFLEX 
flux limit. (This figure is available in colour in the online version of 
the journal.) 

The parent sample is the REFLEX catalogue 
(Bohringer et al., 2004). To ensure the best quality for 
potential targets, the REFLEX catalogue was first screened to 
include only objects which had (i) a firm detection threshold 
of more than 30 source photons in the ROSAT All Sky Survey 
and (ii) a low column density («// < 6 x 10^'^ cm"^). 

Since the Legacy Project selection was intended to be rep- 
resentative of an X-ray flux- or Lx-limited sample, clusters 
were chosen purely on the basis of X-ray luminosity. Further 
selection criteria included: (i) redshift z < 0.2 to sample the 
nearby Universe; (ii) close to homogeneous coverage of the 
luminosity space; (iii) a flux limit corresponding to kr > 2 
keV, to sample the mass range from poor systems to rich clus- 
ters; (iv) detectable with XMM-Newton to approximately a ra- 
dius of Rsm, with distances selected to optimise R500 in the 
XMM-Newton field of view. 

To best assess the scaling relations, the sample should 
have close to homogeneous coverage of luminosity space. The 
luminosity-redshift space was thus sampled in eight almost- 
equal luminosity bins (see Fig. 1)^. The lower redshift bound- 
ary of each bin was placed above the flux hmit curve or close to 
the curve defining the redshift at which Rsoo corresponds to 9' 
(10' for the most luminous clusters). The upper redshift bound- 
aries were defined by the number of clusters to be included in 
the bin (4 objects). For the lowest luminosity bin these criteria 

' One extra bin, containing the most luminous cluster, uses data 
from the XMM-Newton archive. 
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were relaxed, and the bin put at lower redshift, because these 
clusters are fainter. The distribution of clusters in luminosity- 
redshift space is shown together with the luminosity bins cho- 
sen for the Legacy Project sample in Fig. 1. 

Observations of the full sample of 3 1 clusters plus 2 archive 
observations have now been completed. As detailed below, the 
quality of 18 of these observations is not sufficient to derive 
accurate radial temperature profiles to relatively large radius. 
The remaining 15 clusters with good quality data are discussed 
in this paper. As shown in Fig. 1, only one redshift bin is not 
represented in this subsample, thus it is representative of the 
sample as a whole. 

3. XMM data analysis 

Observation data files (ODFs) were retrieved from the 
XMM archive and reprocessed with the XMM-Newton Science 
Analysis System (SAS) v6.1 using the publicly-available cal- 
ibration current as of February 2005. The resulting calibrated 
EMOS and EPN event files were then used in all subsequent 
analysis. 

3.1. ODF preparation 

The data were cleaned for periods of high background due to 
soft proton solar flares using a two stage filtering process. A 
light curve was first extracted in 100 second bins in the [10- 
12]/[12-14] (EMOS/EPN) energy band. A Poisson distribution 
was fitted to a histogram of this light curve, and +3o" thresh- 
olds calculated. A Good Time Interval (GTI) file was produced 
using the upper threshold, and the event list was filtered ac- 
cordingly. Since i) flares often appear to have soft 'wings'; ii) 
the statistics at high energy are often poor; and iii) softer flares 
exist, the event list was then re-filtered in a second pass. In this 
case light curves were made in 10 sec bins in the full [0.3-10] 
band, the smaller bin size being possible because of the greatly 
improved statistics. A histogram was calculated, a Poisson dis- 
tribution fitted, GTIs generated, and event lists were filtered as 
above. 

This type of flare filtering is sufficient in the majority of 
cases. However, it is not as effective in removing flares in cases 
of data sets with softly- varying count rates or gradually increas- 
ing or decreasing count rates. The histograms of such data sets 
invariably have a Poisson distribution with a tail, which is not 
well fitted with a single component. All light curve histogram 
fits were thus carefully examined before further analysis. In 
problem cases, the data sets were cleaned by hand (generally 
by estimating the flare periods by eye, and excluding them). 

Since the object of the present work is to obtain relatively 
high-quality temperature profiles, we only use those observa- 
tions with a cleaned EPN exposure time greater than 10 ks. 
Of the 33 clusters in the Legacy Project sample, 15 meet this 
criterion at the present time."*. These observations are listed in 
Table 1 . Figure 1 shows the X-ray luminosity - redshift distri- 
bution of the LP clusters. Filled circles show the clusters dis- 



Table 1. Basic cluster data. 



RXCJ 


Tx" 




z 




Exp. ^ 


Comments 


0003+0203 


3.71+0.09 


0.085 


4.7 


26,26,17 


A2700 


0020-2542 


5.74 


±0.13 


0.141 


2.2 


16,16,11 


A22 


0547-3152 


6.59 


±0.12 


0.148 


2.1 


23,24,17 


A3364 


0605-3518 


4.68 


±0.11 


0.139 


4.5 


22,23,14 


A3378 


1044-0704 


3.56 


±0.05 


0.134 


3.6 


26,26,18 


A 1084 


1141-1216 


3.60 


±0.08 


0.120 


3.2 


28,28,22 


A 1348 


1302-0230 


3.60 


±0.08 


0.085 


1.7 


25,25,16 


A 1663 


1311-0120 


8.45 


±0.12 


0.183 


1.8 


36,37,29 


A1689 


1516+0005 


4.34 


±0.07 


0.120 


5.4 


26,27,21 


A2050 


1516-0056 


3.75 


±0.10 


0.120 


5.4 


29,30,22 


A2051 


2023-2056 


2.83 


±0.08 


0.056 


5.4 


17,18,10 


S868 


2048-1750 


3.96 


±0.08 


0.148 


4.7 


25,25,19 


A2328 


2129-5048 


3.84 


±0.10 


0.080 


2.2 


21,22,11 


A3771 


2217-3543 


4.60 


±0.08 


0.148 


6.6 


24,24,17 


A3854 


2218-3853 


5.84 


±0.17 


0.141 


5.7 


21,22,11 


A3856 



" Spectral temperature in the radial range 0. 1-0.4 i?200> in keV, es- 
timated using the R-T relation of Arnaud et al. (2005) - see Sect 4.1 
for details. 

* Column density in units of 10"" cm"^ (see text for details). 
^ Cleaned exposure time of EMOSl, EMOS2 and EPN in kilosec- 
onds. 



cussed in this paper, while open circles show clusters for which 
data are pending. REFLEX clusters appearing in the boxes in- 
clude those for which the X-ray criteria (minimum 30 source 
photons and column density «// < 6 x 10^° cm"^) were not met. 
The current subsample is clearly representative of the whole 
sample; only one redshift bin is not represented. 

After removal of periods of high soft proton flux, events 
were filtered according to PATTERN and FLAG criteria. For 
EMOS event files singles, doubles, triples and quadruples were 
selected (PATTERN < 13), while for EPN data sets singles 
and doubles were selected (PATTERN < 5). Events not cor- 
responding to these criteria were removed from the event files 
before further processing. In addition, for all cameras events 
next to CCD edges and next to bad pixels were excluded ( 
FLAG==(5). 

To correct for vignetting, a WEIGHT column was added to 
each event list using the SAS task evigweight. All subsequent 
science products were extracted from this column as described 
in Arnaud et al. (2001). 

Serendipitous and point sources were detected in a broad 
band ([0.3-10.0] keV) coadded EPIC image using the SAS 
wavelet detection task ewavdetect, with a detection thresh- 
old set at 5cr. Detected sources were excluded from the event 
file for all subsequent analysis. 



3.2. Background preparation 
3.2.1 . Flare cleaning 



The remaining poor quality data are being re-observed under 
XMM A04 and A05. 



The basic background files used are those of Read & Ponman 
(2003), which have nominal exposure times of ~ 1 Ms/400 ks 



4 



G.W. Pratt et al.: X-ray cluster temperature profiles 



(EMOS/EPN)^. Close inspection of the high-energy band Hght 
curves showed that considerable periods of high soft proton 
flux still existed in the event files. These periods were removed 
by two applications of the double pass filtering procedure de- 
scribed in Sect. 3.1 above. The resulting filtered background 
event lists have light curve histograms which are adequately 
described by a standard Poisson distribution. The loss in ex- 
posure time is ~ 200/100 ks (EMOS/EPN); the larger relative 
EPN time loss reflects the greater EPN sensitivity to flares. 

After flare filtering, the same PATTERN and FLAG selections 
as above were applied to the background event files. The back- 
ground files were then corrected for vignetting via the addition 
of a WEIGHT column to the event lists. 

3.2.2. Exposure correction 

Since the background data sets consist of stacked observations 
with sources removed, exposure times can vary by up to a fac- 
tor of two across the detector. Using the exposure maps sup- 
plied by Andy Read a new exposure map was computed for 
each event list taking into account exposure variations due to 
the point source subtraction. These exposure maps were renor- 
malised to the new exposure time of the background files after 
flare cleaning. An EXPOSURE column, containing the expo- 
sure time at the position of the event, was then added to each 
background file. The WEIGHT column was then corrected for 
exposure variations by simply dividing by the EXPOSURE col- 
umn. 

3.3. Background subtraction 

The blank sky backgrounds are recast onto the sky using the 
aspect information from the cluster pointing, enabling extrac- 
tion of source and background spectra from the same detector 
regions. This procedure is necessary because the spatial distri- 
bution of the various instrumental Unes is not constant across 
the field of view. 

3.3.1 . Quiescent background 

The XMM-Newton EMOS and EPN backgrounds are domi- 
nated by charged-particle events above ~ 2 keV. The intensity 
of this component can vary by typically +10%, and must be ac- 
counted for by renormalisation. The renormalisation factor for 
each observation was calculated in the source-free [10-12]/[12- 
14] kev (EMOS/EPN) energy band, and flie WEIGHT column 
of each background file was adjusted accordingly. This renor- 
malisation assumes that the particle induced background can 
simply be scaled depending on the count rate. The renormali- 
sation factors are listed in Table 2. 



' The EPN event list is in Extended Full Frame mode and does not 
necessarily contain the same fields as the corresponding EMOS event 
list, hence the shorter exposure time. 

* ftp : //ftp . sr .bham. ac . uk/pub/xituti/expmap* . fits . gz 



cu 




Channel energy (keV) 



Fig. 2. The residual spectrum of R0547, indicating oversubtraction of 
the soft X-ray background (see text for details). Black: EMOSl; red: 
EM0S2; green: EPN. The fit is a Solar abundance MEKAL model 
with kr = 0.24 keV and negative normalisation. The EPN model has 
an additional power-law component with positive normalisation. {This 
figure is available in colour in the online version of the journal. ) 

3.3.2. Soft diffuse X-ray background 

The observation and blank field event files contain a component 
due to the soft diffuse X-ray background, which dominates the 
flux below ~ 1 keV. This component is variable across the sky, 
and thus from pointing to pointing. Correction for this variation 
is thus needed. 

As discussed above in Sect 2, the present cluster sample 
was explicitly defined so that a certain fraction of the detector 
area is essentially free of cluster emission, enabling a direct es- 
timation of the local background in the outer regions of the field 
of view. For each observation, we extracted surface brightness 
profiles in the [0.3-3.0] keV band for each camera. EMOS and 
EPN surface brightness profiles were background subtracted, 
coadded, and binned to 3cr significance. The local background 
spectrum was built using all events outside the radius at which 
the surface brightness profile in no longer significantly de- 
tected. A renormaUsed spectrum from the same region of the 
blank sky background was then subtracted, yielding a residual 
spectrum. We fitted the residual spectrum in the [0.5-10.0] keV 
band with an unabsorbed, solar abundance MEKAL model. 
Energies with significant instrumental line emission (1.4-1.6 
keV for all instruments; 7.45-9. keV for EPN) were excluded 
from the fits. In case of a significant excess of counts in the 
> 2 keV band, presumably due to a remaining component of 
soft protons, a power-law was added to the model. The EPIC 
spectra were fitted simultaneously, with the temperature of the 
MEKAL model linked between instruments. The power-law 
slope and normalisations of all components were free to vary 
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in the fitting process. The normalisation of the MEKAL model 
was allowed to be negative, to account for over-subtraction. 
This purely phenomenological model is capable of describing a 
wide variety of residual spectra (Fig. 2), although care must be 
exercised in deriving the initial model parameters. This residual 
model, with all parameters fixed and the normalisation scaled 
appropriate to the ratio of extraction region areas, was treated 
as an additional component in all subsequent annular fits. 

3.4. Spectral fitting 

If the cluster exhibited an obvious bright cooling core region, 
spectra were accumulated in annuli centred on this surface 
brightness peak. Some clusters have no obvious central peak; 
in these objects spectra were centred on the emission centroid 
evaluated in a 6 arcminute radius. Annular regions for spectral 
fitting were then defined i) to have 1500-2500 EMOSl counts 
available after background subtraction, and ii) to have a mini- 
mum width of 30" to minimise PSF effects. The spectra were 
extracted using the WEIGHT column, assuring full vignetting 
correction. Effective area and response files corresponding to 
the on-axis position were generated using arfgen and rmf gen, 
respectively. The spectra were binned to 3cr significance after 
background subtraction, to allow the use of Gaussian statistics. 
Spectral fits were undertaken in the 0.5-10 keV energy range, 
excluding the 1.4-1.6 keV band (due to the Al line in all three 
detectors), and, in the EPN, the 7.45-9.0 keV band (due to the 
strong Cu line complex). 

Spectra were fitted with absorbed MEKAL models with 
abundances from the data of Anders & Grevesse (1989). The 
residual model described above, with all parameters fixed and 
the normalisation scaled appropriate to the ratio of extraction 
region areas, was treated as an additional component. After 
first checking whether the X-ray absorption was in agreement 
with the HI value, the absorption was fixed at either the HI 
value or the best-fitting X-ray value. The EPIC spectra were 
fitted simultaneously, with temperatures and metallicities tied 
and the EPN spectral normalisation as an additional free pa- 
rameter. Annuli with abundance uncertainties dZ/Z > 0.3 were 
frozen at the average value of the two preceding fitted annuli. 
This procedure will not affect the temperature estimates in view 
of the generally flat abundance profiles in the outer cluster re- 
gions (e.g., De Grandi & Molendi, 2002). 

To take into account systematic uncertainties, the spec- 
tra were initially fitted with nominal background normalisa- 
tion. They were then re-fitted with the background normali- 
sation fixed at +10% of nominal. The changes in the best fit- 
ting cluster temperature were treated as the corresponding sys- 
tematic uncertainty; these were added in quadrature to the sta- 
tistical uncertainties of each annulus. Since the temperature 
determination is dominated by the exponential cutoff of the 
Bremsstrahlung slope at higher energies, which will depend 
strongly on the scaling of the particle background, we believe 
that this approach is extremely conservative in terms of error 
determination. 

Figure 3 shows the observed background subtracted spec- 
trum of the outermost annulus of RXC J2048-1750 (6f75 < 




Channel energy (keV) 



Fig. 3. Observed (background subtracted) spectrum of the outermost 
annulus of RXC J2048 -1750 (6:75 < R < 8:9) The solid line is the 
best fitting model spectrum consisting of a cluster component plus a 
Galactic component (see text for details). (This figure is available in 
colour in the online version of the journal. ) 



R < 8^9). The signal to noise of this spectrum is typical of that 
in the outer annulus across the sample. The solid line shows 
the best fitting model spectrum consisting of a cluster com- 
ponent plus a Galactic component. The fit is excellent, with a 
xl — 0.98 for 157 degrees of freedom. 

3.5. X-ray images 

We produced images for each cluster to enable readers to judge 
the morphology of the 15 objects in the sample. Images of the 
source and associated background files were extracted from the 
WEIGHT column of the EMOS event files in 3'.'3 bins in the 
[0.5-2.0] keV band. (We do not use the EPN for image gener- 
ation due to severe problems with artifacts caused by the large 
gaps between CCD chips in this detector) EMOS 1 and EMOS2 
images were exposure corrected and background subtracted 
separately, after which they were coadded. The total EMOS 
image was then binned to Scr significance using the weighted 
Voronoi tesselation method of Diehl & Statler (2005). 

4. Results 

Cluster images and projected temperature and abundance pro- 
files are described in detail in Appendix A. It is clear that there 
is a general trend for the cluster temperature profiles to decline 
with distance from the centre. For a better understanding of 
how similar, or otherwise, the profiles are, it is instructive to 
look at the scaled temperature profiles. 
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Table 2. Columns: (1) Cluster name; (2) Radius beyond which external region spectra were accumulated (EM0S1,EM0S2,EPN); (3-5) 
Normalisation factor for background rescaling. This factor was calculated from the ratio of the count rate in the observation and background 
files in the [10-12]/[12-14] keV (EMOS/EPN) energy band; (6) Temperature of MeKaL model used to describe the residual spectrum; (7-9) 
xsPEC normalisation of the MeKaL model used to describe the residual spectrum, in units of 10""'. 



RXCJ 



^ext 
(') 



Normexi 



kre, 



MeKaL norm 



EMOSl EMOS2 EPN 



EMOSl EMOS2 EPN 



0003 +0203 
0020 -2542 
0547 -3152 
0605-3518 
1044 -0704 
1141 -1216 
1302 -0230 
1311-0120 
1516+0005 
1516-0056 
2023 -2056 
2048-1750 
2129 -5048 

2217- 3543 

2218- 3853 



11,11,11 

11,11,12.5 

11,11,11 

11,11,11 

10,10,10 

11,11,11 

11,11,11 

11,11,11 

12.5,12.5,13.5 

12.5,12.5,12.5 

11,11,11 

13.5,13.5,14 

11,11,12 

11,11,12 

11.5,11.5,12.5 



1.07 
1.00 
0.99 
1.22 
1.14 
1.06 
1.04 
0.97 
1.16 
0.98 
0.98 
0.96 
1.33 
0.93 
1.24 



0.98 
0.98 
0.93 
1.18 
1.11 
1.04 
1.02 
0.93 
1.14 
0.96 
1.17 
0.99 
1.34 
1.14 
1.25 



1.12 
1.02 
1.04 
1.19 
1.16 
1.13 
1.07 
1.05 
1.29 
1.07 
1.16 
1.06 
1.32 
1.16 
1.25 



0.23 
0.10 
0.24 
0.26 
0.26 
0.27 
0.27 
0.19 
0.24 
0.25 
0.23 
0.20 
0.33 
0.65 
0.26 



-1.01 
-14.32 
-1.38 
-1.74 
-2.16 
-1.22 
-0.81 
0.95 
0.48 
0.75 
1.00 
0.17 
-0.51 
-0.75 
-1.32 



-2.98 
-20.02 
-0.31 
-1.09 
-0.67 
-1.26 
-0.97 
0.95 
1.27 
1.04 
1.16 
0.43 
0.90 
-0.45 
0.30 



-0.54 
-6.00 
-0.76 
-0.42 
-5.00 
-0.40 
-0.41 
0.95 
1.01 
1.67 
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Fig. 4. Scaled projected temperature profiles. Left panel: linear x-axis; right panel: logarithmic x-axis. The profiles have been normalised to 
the mean spectroscopic temperature in the 0. 1 i?200 < r < 0.4i?200 region, where /?2oo has been determined iteratively using the R-T relation 
of Arnaud et al. (2005). The shaded grey area corresponds to the region enclosed by the mean plus/minus the Icr standard deviation. The solid 
line in the left-hand panel is the linear fit in the radial range 0.125 < R200 < 0.5 detailed in Eqn. 1. {This figure is available in colour in the 
online version of the journal. ) 



4.1. Scaled temperature profiles 

We normalise the radial temperature profile of each cluster by 
a global temperature, Tx, which should be representative of the 
'virial' temperature of the cluster. Strong cooling core clusters 
have central temperature decrements of up to a factor of three 
which, when combined with the dependence of the X-ray 
emission, means that average integrated temperatures of such 



systems can be biased. However, the cooling core region rarely 
extends beyond ~ 0.1 Rim- In addition, our measured tempera- 
ture profiles do not extend to much further than 1 Mpc even in 
the best cases, which corresponds to ~ 60 per cent of 7? 200 for 
a 5 keV cluster (Arnaud et al., 2005). We thus chose to use the 
overall spectroscopic temperature in the 0. 1 R2Q0 < r < 0.4 /?200 
region. We estimated this region in an iterative fashion, using 
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the R2Q0-T relation of Arnaud et al. (2005) and starting with 
the mean temperature from the measured temperature profiles. 
The measured values of Tx are given in Table 1 . 

The resulting scaled temperature profiles are shown in 
Figure 4. It is obvious that, despite the large variety of ob- 
jects in this sample, from strong cooling core objects to highly 
unrelaxed systems, there is some similarity in the temperature 
profiles; the profiles generally decline from the centre to the 
outer regions. As an initial measure of the scatter in scaled tem- 
perature profiles, we estimated the dispersion at various scaled 
radii in the range 0.0125-0.5 Rooo- The shaded region in Fig. 4 
shows the region enclosed by the mean plus/minus the Icr stan- 
dard deviation. Clearly the scatter increases towards the central 
regions. The relative dispersion in scaled profiles remains ap- 
proximately constant at ~ 10 per cent beyond 0.1 /?20o- In the 
core regions, however, this increases to ~ 25 per cent. Since the 
profiles have not been corrected for PSF and projection effects, 
this figure is likely a lower limit. In fact there is a clear differ- 
ence between the cool core clusters, which have a large tem- 
perature drop toward the centre, and the non-cool core clusters, 
which generally have profiles which increase linearly or flatten 
toward the centre. 

The largest cluster samples were assembled from 
ASCA and BeppoSAX (Markevitch et al. 1998 and 
De Grandi & Molendi 2002, respectively). More recent 
investigations with XMM-Newton (Piffaretti et al., 2005) 
and Chandra (Vikhlinin et al., 2005) have allowed better 
constraints to be put on the form of the profiles of cool core 
clusters out to 0.4-0.5 /?200- In Fig. 5 we show a comparison 
of our temperature profiles with those from ASCA, BeppoSAX 
and Chandra. We use the same average temperature as above 
to normalise the temperatures, but as in these previous works, 
we scale the radial coordinate to 7?igo using the relation given 
in (Evrard et al., 1996). There is good agreement between the 
profiles measured with different instruments. There is a ten- 
dency for our temperature profiles to scatter around the upper 
edge of the envelope of the ASCA results. We note however 
that the same tendency is seen in the Chandra observations 
(see Fig. 16 of VikhHnin et al., 2005). 

It should be noted that the definitions of the global tem- 
perature and/or the virial radius often differ between samples, 
making exact comparison between different results rather dif- 
ficult. We do not compare with the XMM-Newton results of 
Piffaretti et al. (2005) since their results are quoted for Riso 
measured from the data, rather than derived from the relation 
of Evrard et al. (1996). We also note that the normalisation of 
the Piffaretti et al. (2005) profiles is ~ 20 per cent lower com- 
pared to our results and to those from other satellites. It is 
possible that this difference comes from their different defi- 
nition of global temperature (Piffaretti et al. fit the emission- 
weighted bins outside the cooling core with a constant tem- 
perature). However, the general declining trend with radius is 
similar to their results. 

Fitting the radial range 0.125 < R200 < 0.5 with a simple 
linear model we find, 

T/Tx=Ll9-Q.14R/R2oo (1) 
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Fig. 5. Scaled projected temperature profiles compared with 
the average profiles from ASCA (Markevitch et al., 1998, 
grey band), BeppoSAX observations of cooling core clusters 
(De Grandi & Molendi, 2002, green line), and Chandra observa- 
tions of cooling core systems (Vikhlinin et al., 2005, red line). The 
observed profiles have been scaled using RigQ derived from the 
simulations of Evrard et al. (1996). (This figure is available in colour 
in the online version of the journal.) 



for a fit with the BCES estimator. A Unear least squares fit gives 
identical results. 

4.2. Comparison with simulations 

Negative gradients of the temperature profiles on scales 
R > 0.1/? 200 are naturally produced by cosmological hy- 
drodynamical simulations of galaxy clusters, quite indepen- 
dent of the details of the physical processes included (e.g., 
Evrard etal., 1996; Lewis etal., 2000; Token etal., 2002; 
Borgani et al., 2004; Kay et al., 2004). Markevitch et al. (1998) 
and De Grandi & Molendi (2002) compared observed temper- 
ature profiles from ASCA and Beppo-SAX data, respectively, 
to the results from non-radiative simulations by (Evrard et al., 
1996) and found a reasonable agreement in the outer cluster 
regions. Loken et al. (2002) discussed a universal temperature 
profile in their simulated clusters, whose shape agrees well with 
observations outside the core regions. However, a number of 
authors have shown that including radiative cooling in simula- 
tions causes a substantial steepening of temperature profiles in 
the central cluster regions (e.g., Lewis et al., 2000; Valdarnini , 
2003; Tornatore et al., 2003). The resulting temperature pro- 
files are at variance with respect to the observed properties of 
cool core clusters (e.g., Borgani et al., 2004). This points to- 
wards the need to introduce a suitable energy feedback scheme 
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Fig. 6. Scaled projected temperature profiles compared with the projected profiles of all clusters with kT > 2 keV from the simulations 
of Borgani et al. (2004). The mean profile of the observed and simulated profiles are shown by the black and green lines respectively. The 
observed profiles are scaled by the measured spectral temperature in the 0.1 - 0.4/?200 region. The simulated profiles are scaled using the mean 
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temperature used to scale the profiles. See text for details. (This figure is available in colour in the online version of the journal.) 



to regulate gas cooling in the central regions (e.g., Kay et al., 
2004). 

In Figure 6 we show a comparison of the observed pro- 
files with the projected simulated profiles of all clusters with 
kr > 2 keV from Borgani et al. (2004), in which the SPH code 
GADGET-2 (Springel, 2005) was used to simulate a concor- 
dance ACDM model (Qm = 0.3, Qa = 0.7, erg = 0.8,/; = 0.7) 
within a box of 192 /i"' Mpc on a side, using 480^ dark matter 
and an equal number of gas particles. The simulation included 
radiative cooling, star formation and galactic ejecta powered 
by supernova feedback. The observed profiles are scaled using 
the spectral temperature Tx, measured as described above in 
Sect 4. 1. The simulated profiles were scaled using the emission 
weighted global temperature, with a further 8 per cent adjust- 
ment so that the normalisation in the 0.15 < 7?200 < 0.5 region 
is the same as that of the observed profiles (this adjustment is 
necessary because the emission weighted global temperature 
is not the same as the measured spectral global temperature 
Tx)- Three projections are shown for each cluster. The scatter 
in the simulated profiles is noticeable and comes from colder 
subclumps accreting onto the main cluster and shock fronts due 
to supersonic accretion. The simulated profiles decline contin- 
uously from a peak at about 0.05 Rm)- The mean observed and 
simulated profiles are shown as black and green lines, respec- 
tively. 

There is relatively good agreement in the external regions, 
with the simulated profiles reproducing the observed scatter In 
the central regions the peak of the simulated temperature pro- 
files hes at ~ 0.047?i8o, in contrast to the observations, which 



show a less pronounced peak at ~ 0.06 /?i so, a point which is 
particularly evident from the mean profiles. In addition there 
appears to be considerably more dispersion in the observed 
profiles in the central regions. While admittedly our profiles 
are uncorrected for PSF and projection effects, we note that a 
similar difference in peak position (as compared to the simula- 
tions) is also evident when comparing with the Chandra results 
of Vikhlinin et al. (2005). Clearly, a more rigorous comparison 
would require estimating temperatures in the simulated clusters 
and rescaling their profiles in exactly the same manner as the 
data. Nevertheless, we believe that the agreement between sim- 
ulated and observed clusters is quite good and lends support to 
the capability of numerical simulations to describe the global 
thermal structure of the ICM. 

4.3. Dependence on cluster morphology/dynamical 
state 

It is interesting to investigate whether there are correlations 
between the form of the temperature profile and the mor- 
phology or dynamical state of the cluster. To this end, we 
make a preliminary investigation using the power ratio method 
of Buote & Tsai (1995) to characterise the morpho-dynamical 
state of the objects in our sample. 

4.3.1 . Power ratio calculation 

To obtain a quantitative estimate of substructure in the surface 
brightness distribution of the clusters, we applied the analysis 
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Table 3. Cluster power ratios, calculated in an aperture corresponding to R500, estimated from the R-T relation of Arnaud et al. (2005). 
Columns: (1) Cluster name; (2): Rsoo in arcminutes; (3,5,7) power ratios; (4,6,8) Icr errors on power ratios. 
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Fig. 7. Power ratios. The power ratios are computed in an aperture corresponding to ^500 estimated using the R-T relation of Arnaud et al. 
(2005). See text for details. 



method proposed by Buote & Tsai (1995). In this method, the 
projected mass distribution is associated with the X-ray sur- 
face brightness; a multipole expansion of the X-ray image then 
yields a similar expansion of the gravitational potential. The 
multipole analysis provides a measure of the 'power' of each 
multipole component to the X-ray image of the cluster in abso- 
lute units. To obtain a measure that is independent of the clus- 
ter X-ray luminosity, the 'power' terms are scaled by the zeroth 
order (monopole) moment and are consequently called 'power 
ratios'. The method was recently used to compare substructure 
in nearby and distant cluster samples observed with Chandra 
(Jeltema et al., 2005). 

The method is applied within an aperture radius as de- 
scribed in Buote & Tsai (1995). We used a minimization of the 
first (dipole) moment to obtain an independent centering of the 
cluster emission within the aperture. The lowest order power 
ratios which are of interest for our study are PilPo, P3/P0, 
and PaIP{), which correspond to the quadruple, the hexapole, 
and the octopole moments. Due to the nature of the moment 
functions, a large weight is given to the outer parts within the 
aperture, in particular for the higher moments. Thus the results 



depend somewhat on the choice of aperture. We have explored 
this effect with a range of aperture radii, results from which will 
be described in a future paper. For the purposes of this initial 
investigation, we estimate the power ratios within Rsqo, where 
this radius is estimated using the R-T relation of (Arnaud et al., 
2005). 

Systematic effects and uncertainties for each power ratio 
were also taken into account. We created 1000 simulations of 
each cluster where the image pixels were azimuthally random- 
ized around the predetermined cluster centre. The power ratio 
signal measured for these simulated clusters should be solely 
due to shot noise, providing a measure of the accuracy for re- 
jection of the hypothesis that the cluster has no structure. We 
therefore subtract the mean of the simulated power ratio signal 
from the result obtained for the observed cluster and use the 
dispersion of the simulated results as an approximation of the 
error of the power ratio. We defer a more precise treatment to 
a forthcoming paper The results shown in Fig. 7 demonstrate 
that we have sufficient photon statistics to produce robust re- 
sults with small uncertainties, except for those cases where the 
clusters have a highly symmetric appearance. In general the re- 
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region (a measure of the central temperature dip) with power ratio. 



suits for the power ratios correspond very well to the visual im- 
pression of the cluster, where in particular P2/PQ can easily be 
identified with the cluster ellipticity. The parameter Pj/Pa pro- 
vides then the best measure for further substructure and since 
some ellipticity can be related to the quasi-equilibrium state of 
the cluster, the third moment provides often the strongest indi- 
cation for a deviation from a relaxed dynamical state. 

Power ratio values and Icr errors are listed in Table 3. The 
power ratios for the clusters in our sample in general occupy 
a similar range of values to those derived for a nearby clus- 
ter sample by Jeltema et al. (2005). Plots of P2/A) vs P3/P0, 
P2/A) vs P4/P0 and P3/P0 vs P4/P0 are shown in Fig. 7. Of 
particular note is the strong correlation in P4/P0 vs Pj/Po- 

4.3.2. Preliminary comparison with power ratio 

We can check to see if there are any correlations between the 
power ratio value and the shape of the temperature profile. We 
parameterise the outer temperature slope by measuring the ra- 
tio between the temperature at 0.5 R20Q and the temperature at 
0.27?200 {T{Q.5r,^)/T{0.2r,^). These values are calculated by 
spline interpolation and are extrapolated if necessary. In Fig. 8, 



we show a scatter plot of the outer slope parameter and each 
of the power ratios. We then tested for correlations between 
r(0.5fc„„)/r(0.2R,u„) and power ratio in the linear-log plane. 
We calculate the generalised Kendall's r correlation coefficient 
(Isobeetal., 1986), appropriate for censored data, for each 
scatter plot. The probability that a correlation is not present 
is 80, 66 and 88 per cent for r(0.5R,„„)/r(0.2R,„„) vs P2/P0, 
P3/P0 and P4/P0, respectively. Given the wide range of mor- 
phologies present in the sample, the lack of significant correla- 
tion suggests that, in general, the morphology/dynamical state 
does not have a significant impact on the outer slope of the az- 
imuthally averaged temperature profile, at least at the currently- 
available resolution. 

Another interesting question is whether there is a correla- 
tion between the presence of a central temperature dip and the 
power ratio. In Fig. 9 we show a scatter plot of the ratio of 
the central temperature, (the temperature of the first bin in 
the temperature profile) to the mean spectroscopic temperature 
Tx, and each of the power ratios. Calculating the generalised 
Kendall's t correlation coefficient for each scatter plot, we find 
a probability that a correlation is not present of 55, 87 and 9 per 
cent for Tc/Tx vs P2/P0, P3/P0 and P4/P0, respectively. Thus 
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there is evidence for a weak correlation of central temperature 
drop with P4/P0, in the sense that clusters with smaller P4/P0 
are more likely to have a central temperature drop. We do not 
yet have the full sample of clusters from which correlations can 
be derived, nor have the temperature profiles been corrected for 
PSF and projection effects, which will enhance the observed 
central gradients to some degree. 

It should be noted that the spatial resolution of these and 
other recent temperature profiles, particularly at large radius, 
while much improved over results from previous satellites, is 
still the limiting factor when looking for correlations, or com- 
parison between different cluster subsamples. This will also 
have a bearing on comparisons with numerical simulations. 

5. Conclusions 

We have used XMM-Newton observations of 15 clusters drawn 
from a statistically representative, luminosity-selected sample 
to investigate the behaviour of the temperature profiles. The 
clusters range from morphologically relaxed looking objects 
with strong central surface brightness peaks (e.g., RXC J1044- 
0704), to diffuse structures with significant amounts of sur- 
rounding substructure (e.g., RXC J1516 +0056), and constitute 
a representative subsample. We find that, once scaled appro- 
priately, the temperature profiles are similar in the radial range 
from 0.1 /?200 to 0.5/?200> declining steadily from the central re- 
gions to the outer boundary of the measurements with a relative 
dispersion of ~ 10 per cent out to 0.5 /?200- The region interior 
to 0. 1 /?200 is the region of greatest scatter in the scaled pro- 
files: the relative scatter of ~ 25 per cent is likely a lower limit. 
A preliminary comparison with numerical simulations shows 
relatively good agreement outside ~ 0.1/? 200, with all of the 
measured temperature profiles falling within the scatter of the 
simulated profiles. 

Calculating power ratios for the sample, we investigate 
whether there are correlations between the power ratio mea- 
sured in an aperture corresponding to R^qq and the shape of 
the temperature profile. We characterise the temperature pro- 
file shape in two ways: with the ratio T{0.5 R2oo)/T(Q.2R2oq), 
a measure of the outer slope, and with the ratio TdiT), a mea- 
sure of the central temperature drop. There is no obvious cor- 
relation of outer slope with power ratio; neither is there a cor- 
relation of central temperature dip with P2/P0 or P^/Pq. There 
is evidence for a weak correlation of the central temperature 
dip with PaIPq. The analysis thus suggests that the outer slope 
of the temperature profile is not particularly sensitive to the 
morpho-dynamical state, although there may be some correla- 
tion with the existence of a central temperature drop. Further 
investigation with power ratios evaluated in other apertures, for 
the entire sample, should be undertaken before definitive con- 
clusions can be drawm. 

The overall conclusion from this work on a statisti- 
cally representative sample indicates that clusters are a 
relatively regular population, at least outside the cool core 
regions, with the caveat that comparisons between samples 
or with simulations is limited by the available tempera- 
ture profile resolution. The observed similarity in density 
(Neumann & Arnaud, 1999; Croston et al., 2006) and temper- 



ature profiles (Markevitch et al., 1998; De Grandi & Molendi, 
2002; Piffaretti et al., 2005; Vikhlinin et al., 2005, this 
work) indicates both similarity in the underlying gravi- 
tational mass distribution (such as has already been seen 
in the X-ray mass profiles of morphologically relaxed 
clusters, e.g., Pointecouteau, Arnaud & Pratt 2005), and 
similarity in the entropy of the ICM (such as has been 
seen by e.g., Ponman, Sanderson & Finoguenov 2003; 
Pratt, Arnaud & Pointecouteau 2006). In this case a single 
integrated temperature, excluding the core region, should be a 
good proxy for the total mass. The observed regularity thus has 
important implications for the use of clusters as cosmological 
probes. 

In future papers, we will reinvestigate the trends with the 
full sample, make maps of quantities such as temperature, en- 
tropy and pressure, and estimate the mass, baryon fraction and 
entropy in the clusters. A more extensive comparison with nu- 
merical simulations will also be undertaken. 
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Appendix A: Individual cluster profiles 

A.1. RXCJ0003+0203 

RXC J0003 +0203, also known as Abell 2700, has an average 
temperature of kT - 3.8 keV and lies at z = 0.092. The clus- 



ter presents a symmetric X-ray image but does not possess a 
strong central surface brightness peak. After renormalisation 
of the quiescent background, the spectra extracted in the exter- 
nal region (r > 11') can be adequately fitted with a MeKaL 
model at 0.23 keV with negative normalisation. An additional 
power-law is required for the EMOS2 and EPN spectra. 

The temperature and metallicity profiles are shown in 
Fig. A.l. The temperature profile is consistent with being 
isothermal at large radii, although given the ~ 35% uncertain- 
ties in the final bin, a decline cannot be ruled out. The metallic- 
ity profile is highly peaked, declining smoothly from Z ~ 0.52© 
in the central bin to Z ~ 0.252© at large radii. While the in- 
crease in metallicity towards the centre is reminiscent of a cool- 
ing core (De Grandi & Molendi, 2002), the temperature profile 
does not show a significant central decline, at least at the reso- 
lution of these data. 



A.2. RXC J0020 -2542 

Moderately luminous, lying at z - 0.141 and with a tempera- 
ture of kr = 5.7 keV, this cluster is also known as Abell 22. 
The X-ray image is highly disturbed, with a prominent surface 
brightness edge to the N, and an emission extension to the S. 
The overall X-ray emission is aligned approximately in the di- 
rection of the line joining the two brightest cluster galaxies. The 
annular regions were centred on the surface brightness peak, 
which lies on the surface brightness edge and corresponds to 
the position of the BCG. The residual spectrum shows negative 
residuals and is adequately described with a combination of a 
MeKaL model at 0.10 keV with negative normalisation, with 
an additional power-law component for the EMOS2 and EPN. 

The temperature profile (Fig. A. 2) shows a prominent de- 
cline with radius. The abundances are roughly flat but became 
unconstrained at only ~ 400 kpc from the centre and were 
frozen thereafter. 



A.3. RXCJ0547-3152 

Also known as Abell 3364, this luminous cluster lies at z = 
0.148 and has an average temperature of kT - 6.6 keV. The 
X-ray image shows a bright, offset core, with obvious surface 
brightness edges to the NW and SE. After renormalisation, the 
spectrum of the external region shows negative residuals below 
~ 1 keV. These are adequately described with a MeKaL model 
with negative normalisation and a temperature of 0.24 keV. An 
additional power-law component is needed to fully describe the 
EPN spectrum (see Fig. 2). 

The temperature profile declines monotonically, from kr ~ 
7 keV to kr ~ 4.5 keV, from the centre to the largest radii at 
which we can measure the temperature. The metallicity pro- 
file declines from Z ~ QAZq to Z ~ 0.2Zq between < r < 
300 kpc, but then increases once more to the central value at 
r ~ 600 kpc. The metallicity trends may be connected to the 
disturbed nature of the cluster. 
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Fig. A.l. Image (left), temperature (middle) and abundance (right) profiles of RXC J0003 +0203. In tfiis and the following images, the colour 
scale is square root with a maximum of 0.0015 counts per second (enabling easy comparison between images). The angular size of each image 
has been chosen to approximately match the virial radius of the cluster as determined from the average temperature Tx (as defined in Sect. 4. 1) 
and the R-T relation of Arnaud et al. (2005). The horizontal line without eiTors in the abundance profile plots denotes the regions where the 
abundance was frozen. {Figures are available in colour in the online version of the journal.) 
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Fig. A.2. RXC J0020 -2542. 
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AA. RXCJ0605-3518 

Lying at z = 0.14, with an average temperature of kT - 4.7 
keV, this cluster is also known as Abell 3387. It is highly sym- 
metric, presenting a strongly-peaked central surface brightness 
and no visible substructure. The residual spectrum is well de- 
scribed with a MeKaL model with negative normalisation and a 
temperature of 0.26 keV. An additional power-law component 
is necessary for a full description of the EPN data. 

The temperature profile (Fig. A.4) rises from the central re- 
gions to a peak atR ~ 200 kpc, after which there is a gentle de- 
cline. The abundance profile declines smoothly from ~ 0.6 
in the central regions to ~ 0.2 Zq at R > 500 kpc. The general 
behaviour of the temperature and abundance profiles is very 



reminiscent of the Chandra temperature profiles of cool core 
clusters derived by Vikhlinin et al. (2005). 

A.5. RXC J1 044 -0704 

RXC J 1044 -0704, also known as Abell 1048, lies at z = 0.13. 
It has an average temperature of kT = 3.6 keV and although 
slightly elliptical, is another symmetric cluster with a strongly- 
peaked central surface brightness. The residual spectrum is 
negative in the 0.5-1.0 keV band, indicating oversubtraction of 
the background in this band. The residual spectrum can be fit- 
ted with a MeKaL model at 0.26 keV with negative normali- 
sation. An additional power-law component improves the fit to 
the EPN data. 
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Fig. A.3. RXC J0547 -3152. 
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Fig. A.4. RXC J0605 -3518. 
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The temperature and abundance profiles (Fig. A. 5) are once 
again reminiscent of those of other cool core clusters. The tem- 
perature climbs to a peak at R ~ 200 kpc and declines there- 
after, while the abundance profile declines from the central re- 
gions to the outskirts (although in this case the decline is not 
smooth). 

A.6. RXCJ1141 -1216 

This highly symmetric cluster at z - 0.12 exhibits strongly 
peaked central emission. The cluster is also known as Abell 
1348 and has an average temperature of kT - 3.6 keV. The 
residual spectrum shows negative residuals and is well fitted 
with a simple MeKaL model with negative normalisation and 
a temperature of 0.27 keV. 

The temperature and abundance profiles shown in Fig. A. 6 
are very similar to those of the previous two clusters. The cen- 
tral temperature dip is associated with an abundance enhance- 
ment; the temperature peaks around 200 kpc and dechnes gen- 



tly thereafter The abundance declines smoothly from the centre 
to the external regions. 



A.7. RXC J1 302 -0230 

Also known as Abell 1663, this is a symmetric looking clus- 
ter at kr = 3.6 keV lying at z = 0.085. It has a strong cen- 
tral emission peak. The residual spectrum of the external re- 
gion (r > 11') can be characterised with a MeKaL model at 
0.27 keV with negative normalisation. An additional power law 
component improves the fit for the EMOS2 and EPN spectra. 

The temperature and abundance profiles (Fig. A. 7) are very 
characteristic of cooling core clusters. Compared to similar 
clusters in this sample, however, RXC J 1302 -0230 is charac- 
terised by a particularly steep central temperature drop, and a 
similarly steep increase of metallicity towards the central re- 
gions. 
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Fig. A.5. RXC J1044 -0704. 
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Fig.A.6. RXC Jl 141 -1216. 
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Fig. A.7. RXC J 1 302 -0230. 
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A.8. RXCJ1311 -0120 

This extremely symmetric cluster lying atz = 0.183 is the well- 
known lensing cluster Abell 1689. This is the most luminous 
cluster in the present sample, which is reflected by its particu- 
larly high temperature (kT = 8.5 keV). The residual spectrum 
of the external region (r > 11') shows a positive excess which 
is well modelled by a MeKaL model at 0. 19 keV. An additional 
power-law component improves the fit for the EPN detector. 

The temperature and abundance profiles (Fig. A. 8) are 
not characteristic of cooling core clusters, however. While 
there is a central temperature drop, it is not nearly as steep 
as that displayed by other cool core clusters in this sample. 
Furthermore, the abundance profile does not show a central 
peak. The temperature and abundance profiles we have de- 
rived are in good agreement with those derived (from the 
same XMM-Newton data) by Andersson & Madejski (2004). 
Girardi et al. (1997) use galaxy velocity date to describe this 
cluster as a line of sight merger. This may explain why the clus- 
ters is quite symmetric but does not appear to possess a strong 
cooling core. 

A.9. RXCJ1516+0005 

Also known as Abell 2050, this moderate temperature (kT = 
4.6 keV) symmetric looking cluster lying at z = 0.120 does not, 
however, exhibit peaked central emission. The external residual 
spectrum, accumulated from events from beyond 12.5' from 
the cluster centre exhibits an excess of counts at £ < 1 keV 
and is adequately fitted with a MeKaL model at 0.25 keV. An 
additional power law component improves the EPN fit. 

The temperature profile of this cluster (Fig. A. 9) shows no 
sign of cool core emission, declining linearly from the centre 
to the external regions. As expected, the abundance profile is 
consistent with being flat at Z ~ 0.3 Z© out to 500 kpc (the 
maximum radius at which we can measure abundances). 

A.10. RXCJ1516+0056 

A moderate temperature (kT = 3.75 keV) cluster lying at 
z = 0.120, RXCJ1516-H0056 is also known as A2051. The 
X-ray image presents quite a lot of structure, with several pos- 
sible subclumps at the outskirts of the object. These subclumps 
were excluded from the annuli used to determine the temper- 
ature profile. The background subtracted spectrum of the ex- 
ternal region is well fitted with a single MeKaL model at 0.25 
keV, with positive normalisation. 

The temperature and abundance profiles are shown in 
Fig A. 10. The temperature profile is flat in the inner 400 kpc, 
but declines by ~ 50 per cent at the radius of maximum detec- 
tion. The abundance profile is consistent with being flat out to 
the radius of maximum detection. 

A. 11. RXCJ2023-2056 

Also known as S868, lying at z = 0.056, this is the lowest tem- 
perature cluster in the present sample (kT = 2.7 keV). The ob- 
ject has a fairly regular appearance, but no strong evidence for a 



cooling core. The background subtracted external region spec- 
trum is well described by a MeKal model with positive nor- 
malisation and a temperature of 0.23 keV. An additional power 
law component improves the fit of the EPN data. 

The temperature profile of the cluster, shown in Fig A. 1 1 
is flat in the inner 100 kpc, after which there is a decline. The 
abundance profile declines smoothly in power law fashion from 
a peak of Z = 0.5Zo in the centre to about half that value at the 
outskirts. 

A.I 2. RXCJ2048-1750 

With a temperature of kr = 4 keV and lying at z = 0.085, 
RXC J2048 -1750 presents a fairly disturbed appearance. The 
background subtracted spectrum of the region external to the 
cluster emission can be fitted with a simple thermal model at 
0.20 keV, with positive normalisation. 

The temperature profile of the cluster (Fig. A. 12) declines 
linearly, by more than a factor of two, from the centre to the 
external regions. The abundance profile is very poorly con- 
strained, and we can only measure the three inner bins. 

A.I 3. RXC J21 29 -5048 

A moderate temperature (kr = 3.8 keV) cluster also known as 
A3771, RXC J2 129 -5048 lies atz = 0.08. The X-ray image is 
disturbed, with a distinct elongation in emission from the centre 
towards the NE. The background subtracted spectrum of the 
external region can be fitted with a thermal model at kr = 0.33 
keV with an additional power law improving the fit in all three 
cameras. 

The temperature profile of the cluster (Fig. A. 13) is rela- 
tively flat in the inner 200 kpc or so, but declines beyond this. 
The abundance profile is relatively poorly constrained, but is 
consistent with being flat at an average of Z ~ O.35Z0. 

A.I 4. RXC J2217 -3543 

One of the more distant clusters in the sample, having an av- 
erage temperature of kr = 4.6 keV and lying at z = 0.148, 
RXC J2217 -3543 is also known as A3584. The X-ray image 
is quite compact and symmetric, although the cluster does not 
present strongly-peaked core emission. The spectrum of the 
background subtracted external region presents strongly neg- 
ative residuals below 1 keV and can be fitted with a thermal 
model at 0.65 keV, with negative normalisation. 

The temperature profile shown in Fig. A. 14 declines lin- 
early from the peak of 5.5 keV at the centre to 3 keV at the 
maximum radius at which we can measure the temperature. 
The abundance profile does not show any trends with radius. 

A.I 5. RXC J2218 -3853 

RXC J22 1 8 -3853 is also known as A3856, has an average tem- 
perature of kr = 5.8 keV and lies at z = 0.09. The X-ray image 
is elliptical, presenting an elongation in the SE-NW direction. 
The background subtracted spectrum of the external region is 
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Fig.A.8. RXCJ1311-0120. 
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Fig. A.9. RXC J 1 5 1 6 +0005 . 
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Fig. A.IO. RXC J15 16 -0056 



18 



G.W. Pratt et al.: X-ray cluster temperature profiles 



Arcmin 

4 



Arcmin 
4 



R2023_2056 



R2023_2056 



a7- 



2Dh24mDDs 20h23m3t3B 2Dh23mD0s 20h22m3t3B 2Dh22mDDs 



□ ODDS □ DDI 

(courta'seo'pixel) 



Fig. A.ll. RXC J2023 -2056. 
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Fig. A.12. RXC J2048 - 1750. 




Fig.A.13. RXCJ2129 -5048. 
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well fitted with a thermal model at 0.26 keV, with an additional 
power law component improving the fit for all three cameras. 

The temperature profile (Fig. A. 15) is flat in the inner re- 
gions, rises to a peak at ~ 400 kpc, and then declines (although 
not significantly). The abundance profile is consistent with be- 
ing flat at an average ofZ- Q.3Zq out to 400 kpc, the detection 
limit. 
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Fig. A.15. RXC J22 18 -3853. 



